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We evaluate the crosstalk between adjacent cores in an optical fiber that 
results from electromagnetic coupling. Means of reducing it are discussed. 
We find that a 0.5-nm-thick layer of silver can, in principle, reduce the 
crosstalk from — 20 to — 130 dB without significant increase of the loss. 
These theoretical results are obtained for two identical single-mode dielectric 
slabs. In reality, the slabs are not rigorously identical. Longitudinal fluc- 
tuations of slab thickness reduce the crosstalk by at least 40 dB. The slab 
spacing can accordingly be reduced from, typically, 11 to 6 \xm for a 
constant crosstalk. If the slabs are made dissimilar with a relative differ- 
ence in thickness of 10 percent, the spacing can be reduced further, to 
approximately 1.5 times the slab thickness. For example, a 15-nm spacing 
is required between single-mode dissimilar slabs if the nominal slab thick- 
ness is 10 nm, provided scattering can be neglected. 

I. INTRODUCTION 

In multichannel communication systems, crosstalk between chan- 
nels is a problem that must be considered. Typically, the crosstalk 
should be less than —20 dB. This means that, if an optical power of 
1 mW is fed into one optical guide of a cable, no more than 10 nW 
should be transferred into the other guides. Let us assume a typical 
link length of 10 km. The crosstalk measured over a 1-km-long fiber 
should be less than —40 dB if the power transfer is proportional to 
the square of the fiber length, less than —30 dB if the power transfer 
is proportional to the fiber length, and less than —20 dB if the power 
transfer is independent of the fiber length. As we shall see, the first 
power law is applicable to identical uniform fibers, the second to 
nominally identical irregular fibers, and the third to uniform dissimilar 
fibers. 

In optical fibers, the field decays exponentially in the cladding. 
Therefore, a modest increase in spacing between adjacent fibers is 
usually sufficient to reduce the optical coupling to tolerable values. 
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Yet, in some cases, one needs to minimize the cross section of the cable 
and the spacing between adjacent fibers. Let us briefly discuss a few 
relevant applications. The need for minimizing the distance between 
single-mode cores in a fiber does not arise in communication systems 
presently envisioned for the following reasons: The fiber diameter is 
required to be large (e.g., larger than about 50 j*m) so that the fiber 
is able to sustain mechanical tensions. Thus, quite a few cores can be 
accommodated within the fiber diameter with sufficient spacing. 
Furthermore, the capacity of single-mode fibers is so large there is 
little incentive to introduce more than one core in the same cladding. 
The problem of coupling between single-mode fibers (or between fibers 
carrying few modes) does arise, however, when one tries to increase the 
image-transmission capacity of a fiber bundle up to the diffraction 
limit, each core carrying one bit of image information. Crosstalk 
(image blurring) is minimized if adjacent cores are made dissimilar. 
However, geometrical irregularities may restore a large coupling be- 
tween closely spaced cores. (This, incidentally, raises the possibility 
that measurement of the coupling between dissimilar, closely spaced 
cores gives useful information on the spectral density of the core 
irregularities.) The problem of coupling between single-mode dielectric 
waveguides also arises in integrated optics and in biology in the study 
of the optical behavior of the retina. The results that we present are 
general. They are therefore applicable, in principle, to multimode, as 
well as to single-mode, fibers. However, in practical multimode fibers, 
slow longitudinal variations of the core dimensions make the propaga- 
tion constants of the modes of one core sweep randomly through the 
propagation constants of the modes of the other core. Thus, an aver- 
aging takes place that cannot be ignored. The problem of coupling 
between highly multimoded cores will be only briefly discussed. 

The shielding method discussed in this paper consists of the intro- 
duction of a layer of metal, typically silver, between the adjacent 
optical waveguides. A reservation is in order : In some communication 
systems, metallic layers may be undesirable because they detract from 
the all-dielectric-cable properties. Shielding between adjacent fibers 
can be provided alternatively by low-refractive-index plastics such as 
Teflon® fep (n tt 1.32) that cause the optical field to decay faster 
than in the cladding material. The reduction in coupling, however, is 
much smaller than that provided by metals. Plastic materials can be 
made very lossy by impregnating them with dyes. High losses, how- 
ever, are much less effective than small refractive indices in reducing 
evanescent wave coupling. Therefore, we shall consider mainly 
metallic layers. The practicality of metallic shields remains an open 
question. 
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In the first part of this article series, 1 a general and simple expression 
of the coupling between two lossy open waveguides was derived. Our 
formulation requires that only the normalized fields of the individual 
waveguides along a contour be known. In the present paper, we eval- 
uate in detail the crosstalk between two parallel slabs caused by the 
electromagnetic coupling and means of reducing it. The crosstalk 
between two optical slabs has been evaluated by Marcuse, 2 although, 
in Marcuse's work, the slabs are assumed identical. In reality, un- 
avoidable fluctuations in the slab dimensions reduce the crosstalk, as 
we shall see, by more than 40 dB. Marcuse has also evaluated the reduc- 
tion of crosstalk provided by a layer of absorbing material located 
between the slabs. He found that the waveguide loss increases to in- 
tolerably high values before any significant reduction in coupling can 
be obtained. We find that, if the intermediate layer is metallic, the 
coupling can be drastically reduced without any significant increase 
of the waveguide loss. This discrepancy results from the fact that, for 
metallic layers, the permittivity is negative. For very dissimilar media, 
the first-order perturbation used by Marcuse is not applicable. In the 
present paper, we assume that the perturbation caused by the inter- 
mediate layer on the propagation is small, but we do not assume that 
the field in that intermediate layer is close to the field that would 
exist in the absence of the layer. 

In Section II, we evaluate the crosstalk between optical waveguides 
when the axial wave numbers (or propagation constants) of the iso- 
lated guides fluctuate along the system axis. In Section III, we eval- 
uate the spacing between slabs corresponding to a given crosstalk. In 
Section IV, the transmission is evaluated of a metallic layer under 
evanescent wave excitation and the crosstalk reduction. In Section V, 
we evaluate the loss that results from the introduction of a metallic 
layer near a slab waveguide. In Section VI, a simple approximate 
formula is given for the coupling between oversized round fibers. It 
is compared to exact results. Finally, brief comments are made in 
Section VI concerning the applicability of quasi-ray optics techniques 
in evaluating the coupling between irregular oversized fibers and the 
effect of bending. A few general results that do not seem available in 
convenient form in the literature are derived in the appendices. 

II. FAST COUPLING 

Solution of the coupled-mode equations when the axial wave num- 
bers of the isolated guides are constant, or vary linearly with z, is 
recalled in Appendix A. In the present section, only the results are 
given. 
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Let us first assume that the coupling c between the two guides and 
the axial wave numbers k h k 2 , of the isolated guides is constant 
(independent of z). Let a power unity be fed into guide 1 at z = and 
the other guide, guide 2, be unexcited. The power in guide 2 grows, at 
first, according to the law (see Appendix A) 

P 2 (z) = {cz)\ (1) 

This result is valid only as long as Az « 1, where we have defined 

A a {[(ft, - k 2 ) 2 /4T\ + c 2 }*. (2) 

For example, a -20-dB crosstalk (P 2 = 0.01) over a 1-km length 
of cable is obtained, according to (1), if c = 10~ 4 nr 1 . Condition 
Az « 1 is, for identical guides, z « 10 km. However, if (ki — k 2 )/ 
(ki + k 2 ) = 10~ 4 , law (1) is applicable only if z « 1 mm, a dras- 
tically different condition. In Section III, the distance between the 
guides that corresponds to this particular coupling is evaluated. 

Now let ki — k 2 vary linearly with z. The coupling c remains a 
constant. We write 

ki{z) = k + az, k 2 (z) = k — az, (3) 

where k and a denote constants. At large \z\, the coupling is insig- 
nificant because of the large value of ki - k 2 . The coupling becomes 
important only near the origin, z = 0, where near-synchronism is 
achieved. Let a power unity be fed into guide 1, at large negative z. 
The power transferred to guide 2 at large positive z is exactly (see 
Appendix A) 

P 2 = 1 - exp (-7TC 2 /a). (4) 

We are interested in the case where the A;'s are crossing very rapidly. 
Thus, let us assume that a is large and that, consequently, -r&la is 
small. In that approximation, 

p 2 = wc 2 /a « 1. (5) 

In most practical systems, ki — k 2 oscillates as a function of z. A sig- 
nificant amount of coupling between two guides takes place only near 
the crossing points. To develop an understanding of the effects of 
longitudinal variations of the difference of the axial wave numbers 
ki(z) and k 2 (z), we model the difference in wave numbers as a simple 
sinusoid, i.e., 

i(ki - k 2 ) = S sin (flz), (6) 

where 8 denotes the peak deviation of {ki — k 2 )/2 and 27r/fi the period 
of oscillation. It seems reasonable to assume that the phases of the 
signals picked up by fiber 2 at the successive crossing points are un- 
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correlated and that, consequently, the powers add up. This inco- 
herency is a consequence of the fluctuations of the phase of the optical 
field between successive crossing points. According to (6), the slope 
a introduced in (3) is 

a = fifi. (7) 

The number of crossing points over a length z is Slz/ir. Thus, the power 
collected by guide 2 over length z is 

P 2 = (G*/7r)(7rc 2 )/(5n) = c 2 z/8. (8) 

Note that P 2 is independent of SI. P 2 is proportional to c 2 , as was the 
case in the absence of fluctuations, but it varies linearly with z rather 
than being proportional to z 2 . Let us compare P 2 in (8) and P 2 in (1). 
The ratio of these two collected powers is 

P 2 (uniform fibers) __ . . 

— oz. (v) 



P 2 (nonuniform fibers) 

It seems reasonable to assume that, over a length of 1 km (z = 10 9 
/xm), the relative variations of the axial wave number are larger than 
10 -4 : 8/k > 10 -4 . For the single-mode slab considered in the next 
section, this number corresponds to a fluctuation of the slab thickness 
of 0.01 jum. Because k is of the order of 2-k /xm _1 , the reduction in cou- 
pling owing to the lack of identity between the two slabs is, in that case, 
of the order of 50 dB. The results obtained are therefore much too con- 
servative if we assume that the optical guides are identical in evaluat- 
ing the crosstalk. 

III. EVALUATION OF COUPLING BETWEEN TWO SLABS 

Let us consider two identical dielectric slabs having thickness 2d 
and material free wave number k. The free wave number in the medium 
between the slabs (cladding) is denoted k s , and the spacing between the 
slabs is denoted 2D. (See Fig. 1. The intermediate layer is to be ignored 
for the moment.) The expression for the coupling c between the fun- 
damental H waves is well known (see, for example, Ref. 1) : 

c = xRexp (-2kD), (10) 

where 

k = {k\ - *»)* (11a) 

R = (M)- l [l + (I/**)] -1 !! - (K 2 d 2 /P 2 )] (lib) 

F 2 = (k 2 - k 2 s )d 2 . (lie) 

k z denotes the axial wave number of the isolated slabs (previously 
denoted k\ and & 2 for the two waveguides). If we require that only one 
H mode propagate (for simplicity, we shall ignore the E waves), the 
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Fig. 1 — Coupled dielectric slabs with thickness 2d and free wave number k. The 
cladding medium has free wave number k.. Crosstalk can be reduced by introducing 
a metallic layer with free wave number k m (almost purely imaginary) and thickness 
2d m . 

maximum value of F is ir/2. The theory of dielectric slabs shows that, 
for that value of F, xd = 1.28. Thus, the coupling is 

c = (0.24 /fcd 2 ) exp (-2.5QD/d), (12) 

where we have made the approximation k z tt k in the first term. Thus, 
for a constant relative spacing D/d, the coupling between two single- 
mode slabs varies as the inverse of the square of their thickness. 

Let us evaluate c for the numerical values 

2d = 1.32 /im, k = 2t X 1.45 par 1 , k s = 2t X lApanr 1 . (13) 

Thus, 

F = w/2, k z d ttkd = 8.88. (14) 

If we substitute these results in (lib) and (12), we obtain 

R = 0.021, c(in m" 1 ) = 4 X 10 4 exp (-3.88D), (15) 

where D is in /xm. If the slabs are identical, — 20-dB crosstalk in 1 km 
is obtained, as we have seen in Section II, when c = 10 -4 m -1 . This 
corresponds, according to (15), to a spacing 



2D = 11 Mm. 



(16) 



If the slabs have some irregularities, with 8/k = 10 -4 (corresponding 
to a variation of slab thickness of 0.01 fiva), —20-dB crosstalk is ob- 
tained when c = 0.25 m -1 . This coupling corresponds to a smaller 
spacing : 2D = 6.2 nm. If the slab thickness is chosen equal to 10.5 
/xm, keeping F = v/2 (An/n = 5 X 10 -4 ), the spacing required for 
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identical slabs and — 20-dB crosstalk over a 1-km length is 2D = 66 
/xm, a rather large spacing. 

If the two slabs are uniform but are made deliberately dissimilar, 
a lower crosstalk is obtained. The relative difference 8/k in axial wave 
numbers is approximately 0.5 (Ad/d)/(kd) 2 , where Ad/d is the relative 
difference in thickness of the two slabs (F = t/2). For example, if 
one slab has a thickness 2d equal to 1.32 /xm and the other has a thick- 
ness equal to 1.2 jxm, the relative difference in k z is : 8/k = 0.65 X 10 -3 . 
The maximum relative power that can be transferred from one slab 
to the other is, according to eq. (39), equal to (c/A) 2 , where A ~ 5. 
Thus, a —20-dB crosstalk corresponds, for the above value of 8, to a 
coupling c = 5/10 = 580 m _1 . The slab spacing 2D corresponding to 
that coupling is given by (15). We obtain 2D = 2.2 /xm. More generally, 
we find that D Fzi 1.5d for any d, if F is kept equal to 7r/2 and Ad/d 
= 0.1. Thus, a considerable reduction in spacing is tolerable, in prin- 
ciple, if the slabs are made dissimilar. However, fast fluctuations 
along the z-axis of the slab dimensions with a period of the order of 
ir/8 Pa 100 nm would reestablish synchronism between the two slabs. 
Fluctuations that are too small in amplitude to deteriorate the prop- 
agation under normal conditions (e.g., no significant coupling to the 
radiation modes) may nevertheless introduce a large crosstalk when 
the slabs are very close to each other. Thus, the result obtained above, 
that the spacing between slabs can be reduced to 1.5 X (2d) if the 
fibers are made dissimilar, may not hold true in practical conditions. 

IV. TRANSMISSION THROUGH A METALLIC LAYER UNDER 
EVANESCENT WAVE EXCITATION 

The results in Section II show that the crosstalk power P 2 is pro- 
portional to the square of the coupling c. We have shown in Ref. 1 
that, for identical slabs and a symmetrical configuration, the coupling 
c is proportional to the square of the normalized field halfway between 
the two slabs. Thus,. the crosstalk is proportional to the fourth power 
of the normalized field halfway between the two slabs. If we introduce 
a metallic layer of thickness 2d m , symmetrically centered between the 
two slabs as shown in Fig. 1, the crosstalk is reduced in proportion to 
the fourth power of the field in the middle of the metallic layer. This 
field reduction, denoted t (for transmission), is evaluated in the present 
section. 

Let us consider an evanescent wave with axial wave number k z > k„ 
where k e denotes the free wave number in the medium. This wave 
decays in the x direction according to 

E{x) = E exp (~kx), (17) 

Ktt{k*-k^. (18) 
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Let us now introduce a metallic layer with complex wave number 
k m = k mr + ikmi and thickness 2d m . The ratio t of the field in the 
middle of the layer to the field at the same point in the absence of the 
layer is derived in Appendix B. Provided the layer is sufficiently thick 
or, more precisely, that 



where 
we have 



Real (K m d m ) » 1, (19) 

K m m (*■ - k 2 m )i, (20) 

t = [AkkJ{k + K m )*~] exp [(« - K m )d m J (21) 



At a free-space wavelength X = 1 /mi, A; = 2ir /mi l , the wave number 
of silver is almost purely imaginary, 3 

k 2 m = (k mr + ik mt y = (0.2k + i5k ) 2 = -985 + 79* (in /mi" 2 ), (22) 

and, for a typical glass, assumed lossless (n, = 1.4), 

k 2 s = nlkl = (1.4fc„) 2 = 77.4 /mi -2 . (23) 

With the value of k 2 = k\ - k\ = 3.76 /xm~ 2 in (14), and k 2 m , k\ in (22) 
and (23), we obtain K m = 32 — 1.3*, and, from (21), a power trans- 
mission 

T = tt* = 0.062 exp (-60d m ), (24) 

where d m is in /xm, provided 

d m » 0.03 jum. (25) 

Because the crosstalk power P 2 is proportional to the square of the 
power transmission T, the introduction of a layer of silver of thickness 
2d m between the two slabs reduces the crosstalk in dB by 

201ogio(!T) = 520d m , (26) 

where d m is in /zm. For example, if the layer thickness is 2d m = 0.5 pm, 
the crosstalk is reduced by 130 dB. This reduction is independent of 
the initial value of the crosstalk, within the approximations made. 
Thus, a 0.5-/im-thick layer of silver is sufficient to ensure a complete 
isolation of adjacent fibers, at a wavelength \ = 1 M m - 

Surface polaritons can be guided near the dielectric (k 2 s > 0) and 
metallic (km < 0) interface. However, the losses of such modes are 
extremety high over a distance of 1 km. The cladding modes are also 
strongly attenuated, and it seems that they can be safely ignored. For 
comparison, let us consider, in place of the metallic layer, a low-index 
plastic material of the Teflon type, with a refractive index n = 1.32. 
We now have km = 69 Mm -2 and n m = 3.47 /mi -1 . We obtain a cross- 
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talk reduction equal to 26d m in dB, where d m is in urn. Thus, a 50-dB 
reduction in crosstalk requires a 4-/xm-thick layer of low-index plastic 
material. 

V. LOSS INTRODUCED BY A METALLIC LAYER 

We are now concerned with the fact that, because the refractive 
index of a metal is not purely imaginary, the presence of the metallic 
layer may increase significantly the loss of the modes guided by the 
fibers. This loss depends critically on the distance between the metallic 
layer and the fibers and, therefore, on the distance between the two 
fibers. The loss suffered by the fiber is influenced by the complex re- 
flection of the metallic layer for evanescent waves. This reflection, 
strictly speaking, depends on the thickness of the metallic layer. Exact 
expressions are given in Appendix B. However, in all our numerical 
examples, the thickness of the metallic layer is so large that it can be 
assumed infinite. In that case, the reflection r reduces to 

r = (k - K m )/{K + Km ), (27) 

where k and n m are denned in (18) and (20), respectively. Because the 
imaginary part *„„• of k„, is much smaller than the real part K mr , the 
imaginary part r, of r is approximately 

r,- « 2/cK rai //cL. (28) 

If we use for k z , k„ and k m the numerical values in (14), (23), and (22), 
respectively, we find r, = 0.005. 

To obtain the loss suffered by the slab, we use the perturbation 
formula derived in Appendix C. The variation of k z is assumed to be 
small. The variation of the field near the perturbing object, however, 
is not assumed small. In the present case, k z is real before perturbation. 
The introduction of the metallic layer causes k, to acquire a small 
imaginary part, k zi . The imaginary part k zi of k z is the fiber loss, in 
neper/unit length. There is also a small variation of the real part of 
k z . This variation, however, is of no interest to us. We have (see 
Appendix C) 

k zi = rudi exp (-2/cD m ), (29) 

where R is the slab parameter defined in (12a) and D m the distance 
between the slab and the metallic layer. The imaginary part n of the 
metallic layer reflectivity is given in (28). 

For the numerical values used earlier in (14) and (15), we obtain 
from (29) 

£d B/km = 8.7 X 10 3 A- ZI = 2.6 X 10 6 X exp (-3.88D m ), (30) 
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Fig. 2— Reduction in crosstalk and increase in fiber loss resulting from the intro- 
duction of a silver layer of thickness 2d m (free-space wavelength = 1 /*m). The dielec- 
tric slabs have a normalized frequency F= (fc 2 - k])>d = x/2. Their spacing is kept 
equal to 11 /*m. The loss varies with d m only because of the change in the slab-layer 
spacing. In the absence of metallic layer, crosstalk is —20 dB/km. 

where D m is in pm. For D m = D - d m = 5.25 Mm, the loss introduced 
by the metallic layer, given in (30), is only 

£ = 0.017 dB/km. (31) 

This loss is quite negligible compared with the other losses (absorption 
because of impurity or scattering losses) suffered by the wave. How- 
ever, because £ depends critically on D m , this loss may not be negligible 
in all practical cases. The reduction of the crosstalk and the increase 
of loss caused by a silver layer of thickness 2d m are shown in Fig. 2 for 
the dielectric slabs considered earlier, as functions of 2d m . Note that, 
if we assume for simplicity that the thickness of the metallic layer is 
negligible compared with the slab spacing (2d m « 2D), the (dimen- 
sionless) ratio of k zi (loss) and c is, approximately, 



k zi /c = 2(k mr /k 2 mt )K*d. 



(32) 



Thus, the best metal, from the point of view of propagation, is the 
one whose k mr lk 2 mi is the smallest. 

VI. ROUND FIBERS 

The general coupling formula in Ref. 1 is applicable, in principle, 
to round fibers. Round fibers are more often encountered in practice 
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than are slabs. The geometry is shown in Fig. 3. The fibers are assumed 
identical, with radius a and spacing 2D. The results are given only for 
the scalar fundamental field ^ £^ HE U of oversized fibers [F = 
(k 2 — k 2 )*a » 1]. In that approximation, the normalized field is 
easily found to be (see Part II of Ref. 1) 



#(y) = UoU'khF)- 1 exp (-Fy 2 /2a 2 ), 



(33) 



where wo ~ 2.4 • • • is the first zero of the Bessel function of order zero. 
The y axis is tangent to the rod considered, as shown in Fig. 3. The 
Fourier transform of if/(y) is 



ip(k y ) = (2t)-*£%(2/) exp (-ik y )dy 
= 7r-W-*F-i exp (-k 2 a 2 /2F). 



(34) 



Because the spectral component \f/(k u ) varies approximately as 
exp (— sx) as a function of x, where s= (k\ — k 2 ) J £^ F/a, the 




Fig. 3 — Variation of the coupling between two dielectric rods of radii a as a func- 
tion of their spacing (2D). Tne dimensions and free wave numbers are shown. The 
parameter K is denned as ca(l — tj/fc 2 ) -1 , and c is the coupling. The plain lines are 
from Snyder exact theory, 7 and the dashed lines from the theory in Ref. 1, applied 
to large normalized frequencies F. 
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coupling is 

S(ky)$ 2 (ky)$i(ky)dky 
-00 

= (u 2 /ir l >)k- 1 F-ia- 2 exp(-2FD/a). (35a) 

In place of c, we can use a normalized coupling K defined by 

K = ca[_l - (^A 2 )]->. (35b) 

In the general expression for c in (34), ^ 2 and $i represent the spectral 
components of the field of the two fibers along the y-axis at x = 0. 
The normalized coupling K is plotted in Fig. 3 (dashed lines) as a 
function of the ratio D/a of the fiber spacing (2D) to fiber diameter 
(2a). In that figure, the parameter is the normalized frequency F. For 
comparison, an exact result obtained by Snyder 4 is shown as a plain 
line. The agreement is very good for F > 4. 

The advantage of the method used in this section is that it is ap- 
plicable when the two fibers are separated by a metallic layer. In that 
case, one need only introduce inside the integral sign in the first ex- 
pression in (34) a term T(k u ), where T denotes the power transmission 
of the metallic layer, defined in (24). T now depends slightly on k y 
because, in the expressions given earlier for T, the axial wave number 
k z should be replaced by (k 2 z + k 2 ,)*. The effect of the dependence of T 
on k u is small, however, and the value obtained earlier for T for slabs 
is approximately applicable to round fibers as well. 

VII. MULTIMODED IRREGULAR FIBERS 

We shall make only qualitative comments. In the preceding calcula- 
tions, we have considered the coupling between one mode of one core 
and one mode of another adjacent core. If the cores can carry many 
modes and have dimensions that fluctuate as a function of z, with such 
an amplitude that the variations in axial wave numbers exceed the 
spacing (in axial wave numbers) between adjacent modes, some averag- 
ing takes place. The situation becomes comparable, at least over some 
distance, to that of a slab radiating power into a semi-infinite dielectric, 
a situation discussed in detail in Part II of this series of papers. 1 

Let us picture the field in slab 1 (excited at z - 0) as made up of 
two plane waves. The plane wave moving toward slab 2 tunnels into 
slab 2. Because of the fluctuations in axial wave numbers, the power 
transferred from slab 1 to slab 2 is essentially the power carried by 
that tunnelling wave; we can ignore the fact that this wave, after 
tunnelling, is reflected back and forth inside slab 2 and may tunnel 
back to slab 1. The power transferred from slab 1 to slab 2, then, is 
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proportional to z, for small z, rather than to the square of z, as is the 
case in the absence of irregularities. This picture is consistent with 
that used by Cherin, 5 who adds the powers transmitted by tunnelling 
rays. Let us emphasize that the validity of this quasi-ray optics ap- 
proach rests on the presence of large slow fluctuations of the core 
dimensions. A simple calculation shows that the relative fluctuations 
of the slab thickness must exceed the reciprocal of the mode number. 
This condition is never met for the low-order modes, but it may be 
met by the higher-order modes. Thus, the situation is rather compli- 
cated and requires a deeper analysis. This quasi-ray technique should 
not be confused with that of Kapany and Burke, 6 where the slabs 
are assumed identical and the fields of the tunnelling rays, rather 
than their powers, are added. In the preceding discussion, we have 
assumed that the fiber cable is essentially straight. The coupling 
increases significantly if the cable is bent. 7 This effect makes it even 
more important to provide shields between adjacent fibers. 

VIII. CONCLUSION 

We have shown that a drastic reduction of crosstalk between parallel 
dielectric slabs can be obtained by introducing a layer of silver (thick- 
ness f^ 0.5 jum) between adjacent slabs. The reduction, in decibels, is 
proportional to the imaginary part of the refractive index of the metallic 
layer and to the layer thickness. In many cases of practical importance, 
the loss introduced by this metallic layer is negligible. We have also 
shown that, because of unavoidable irregularities in the fiber dimen- 
sions, the crosstalk is at least 40 dB below that expected for identical 
fibers. 
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APPENDIX A 

Fast and Adlabatic Coupling 

Let $ denote the field of a guide, such that W* is the power. When 
two guides are weakly coupled, their respective fields \pi, \f/2 approxi- 
mately satisfy the well-known equations 8 

-idfi/dz = &i(z)^! + cip 2 , . 

-idfa/dz = k 2 (z)f 2 + cfa. ( ' 

For simplicity, we assume that the axial wave numbers fc x , k 2 of the 
isolated guides are real and that the coupling c is a real constant. The 
solution when hi, k 2 are constant is well known. For the convenience 
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of the reader, this solution is derived below. The general solution of 
(36) is a superposition of normal modes 

yp\(z) = $t exp (ik+z) + ^f exp {ikrz) 

}p 2 (z) = }p£ exp (ik+z) + fa exp {ikrz), 

where 

fc± = (ki - k 2 )/2 ± A 



(37) 



A = {[(fci-fc 2 )74] + c*}*. 



(38a) 
(38b) 



If the initial conditions are MO) = 1, ^(0) = 0, that is, if only 
guide 1 is excited at the origin (2 = 0), the field in the unexcited guide, 
2, is 

* 2 (z) = (*VA) exp [>'(&i + £ 2 )2/2] sin (Az). (39) 



Thus, for small z, the power in guide 2 increases as 

P 2 (z) = (cz) 2 , Az«l. 
This result is independent of k x — k 2 . See Fig. 4. 



(40) 




Fig. 4— Variation of the optical power picked up by fiber 2, where only fiber 1 is 
excited at z = 0, as a function of the normalized axial distance. The axial wave 
numbers of the isolated fibers are assumed to be constant but different [parameter 
(fci - kj)/2c]. Note that the behavior for small cz is independent of hi — k t . 
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(b) 

Fig- 5 — (a) Linear variation of the axial wave number of the isolated waveguides 
as a function of the axial coordinate z. The hyperbolas represent the normal mode 
wave numbers, (b) Adiabatic coupling in fiber optics. All the power from one fiber 
is transferred to the other fiber if the fc'a vary sufficiently slowly. This principle is 
applicable to multimode fibers. 



Let now the axial-wave numbers k h k 2 of the isolated guides vary 
linearly with z 

ki(z) = k + az, k 2 (z) = k - az, (41) 

where k and a denote constants. Synchronism takes place only near 
the origin, 2 = 0. Let us set 



<Ai.2(z) = Ai,*(z) exp (ikoz) 



(42) 



in (36). After differentiation and substitution, we obtain an equation 
for A i, 

(ffiAi/dz*) + (a 2 * 2 + c 2 - ia)A x = 0. (43a) 

A similar equation holds for A 2 that we need not write down. Equation 
(43a) is the equation for parabolic cylinder functions. The asymptotic 
form of the solution, valid for —r/2 ^ arg(z) ^ *■ is, for a power 
unity at z = - oc (see Ref. 9), 

Ai(g) = exp p(a/2)2 2 + i(c 2 /2a) log (-2)], 2 « c/a (43b) 

Ai(z) = exp [z(a/2)2 2 + i(c»/2a) log (2) - «*/2a], 2 » c/a, (43c) 

as we easily verify by substituting (43b) in (43c) and neglecting terms 
of order z~ 2 . To go from (43b) to (43c), note that log (-2) 
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m fa + i g ( 2 ). Note also that a change in the unit with which z is 
measured affects only the amplitude of Ax, which is arbitrary. 

The power in guide 2 after the interaction has taken place, that is, 
for large positive z, is, according to (43c), 

F 2 - 1 - AiA* = 1 - exp (-7tc 2 /«)- (44) 

Let us first assume that rc 2 /a is very small compared with unity, 
that is, the fc's are crossing very rapidly. In that case, guide 1 transfers 
only a small amount of power to guide 2, equal to tc 2 /<x. This is the 
result used in the text. 

When TC 2 /a is very large compared with unity, that is, when the 
variation of ki - fa is very slow, almost all the power from guide 1 
is coupled to guide 2. This is the principle of the Cook adiabatic 
coupler. 10 This mechanism is applicable also to multimode dielectric 
waveguides. It may be used to couple two optical fibers because the 
dimensions are not critical. Only slowness is required. 11 (See Fig. 5.) 

APPENDIX B 

Transmission and Reflection at a Metallic Layer Under 
Evanescent Wave Excitation 

Let the metallic layer have a complex free wave number k m = k mr 
+ ik mi and a thickness d m . The surrounding medium is assumed to 
have a real free wave number k,. The field has the general form (see 

Fig. 6) 

f£o[exp (— kx) + r exp {kx)~\, 
E(x) = \E~exp (-K m x) + E+exp (k„x) 
l&y exp (— kx) 



x ^ 

^ x ^ d, 

x > d m . 



(45) 



where 



is real, and 



K = {k\ - k*)> 

Km = (jb> - *y*. 



(46) 
(47) 





x 


k. 








K m 









— — — m 



Fig. 6— Transmission of a metallic layer with thickness d m and free wave number 
km under evanescent wave excitation (axial wave number k, > k.). At large negative 
x, the field is assumed unperturbed by the layer. 
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The axial wave number k z is assumed to be real and larger than k„. 
By specifying that E and dE/dx are continuous at the boundaries, 
x = 0, x = d m , we obtain the reflection r and the transmission t: 

r = [(*/«.) - (K m /«)j2coth U m d m ) + {kJk) + («A ffl )]-' (48) 

t = exp («*„){ cosh ( Km d m ) + }[((»/*) + (k/k»)] sinh (kX))" 1 . (49) 

We shall now assume that the metallic layer is thick in the sense 
that Real (ic m d m ) » 1. These conditions are well satisfied for the 
metallic layers that we consider in the main text. In that case, (48) 
and (49) reduce to 

r = (« - «.)/(* + Km ) (50) 

t = [4:kk„,/(k + <c m ) 2 ] exp [(« - K m )d m ], (51) 

respectively. Equations (50) and (51) are the results used in the text. 

APPENDIX C 

Loss Introduced by a Metallic Layer 

Let us consider a uniform reciprocal waveguide and let a uniform 
rod be introduced that perturbs the propagation of the waveguide 
(Fig. 7a). We assume that the perturbing rod does not support trapped 
modes or, if it does, that the axial wave numbers of these trapped 
modes are sufficiently far away from that, k z0 , of the waveguide. No 
resonant coupling is assumed to take place. 

We shall first recall a very general result. Let E + , H + and E p , H p 
denote two time-harmonic fields at the same frequency in the same 
medium. If we assume that the medium is reciprocal (that is, that the 
tensor permittivity is symmetrical), it readily follows from the Maxwell 
equations that the divergence of the vector 

J = E+ X H p + H+ X E p (52) 

is equal to zero. Thus, the flux of J through any closed surface is equal 
to zero. In what follows, an exp ( — iut) term is omitted. 

Now let E+, H + be the field propagating in the —z direction along 
an open waveguide. The dependence of E + and H+ on z is denoted : 
exp ( — ik z0 z). Let E p , H p be the field propagating in the +z direction 
in the presence of the perturbing rod with an exp (ik t z) dependence 
on 2. The closed surface S is taken as the surface shown in Fig. 7a 
bounded by the planes 2 = and z = dz, the volume of the perturbing 
rod being excluded. For that choice, the medium enclosed by S is the 
same for both fields. We can therefore use the result stated earlier that 
the flux of J through S is zero. Let us consider the various contributions 
to that flux. The flux of- J through the plane z = dz differs from the 
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PERTURBING ROD 





(b) 

Fig, 7— ( a ) Schematic for the derivation of the general perturbation formula. The 
dielectric waveguide is perturbed by a small lossy rod. The closed surface S extends 
to infinity where the fields considered are assumed to vanish, (b) Application to the 
perturbation of H waves guided by a dielectric slab (k) by a lossy slab (k m ). 



flux of J through the plane 2 = only by a factor — exp [i(k, - k zQ )dz']. 
The difference between these fluxes is, therefore, i(k, — k z0 )dz for 
small dz. Because we are considering only trapped modes, the flux at 
infinity is zero. The flux through the surface surrounding the perturb- 
ing rod is dz times the line integral of J-dC, with dC a vector perpen- 
dicular to the contour surrounding the rod, pointing inward, whose 
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length is the elementary arc length. Thus, we have exactly 

i&k. = J c J-dc/ J s JdS, (53) 

where 

Ak z m k, - fc, . (54) 

So denotes the transverse plane, 2 = minus the area enclosed by C, 
and dS denotes a vector directed along the z axis whose length is the 
elementary area. The derivation given above is almost identical to 
that in Ref. 1 for coupled waveguides. We now assume that the per- 
turbation is small. Thus, we can replace E p , H p by the unperturbed 
field E, H propagating in the +z direction in the integral over So in 
(53). This is not permissible, however, for the integral over C, in 
general. 

Let (52) be specialized to the H waves guided by a dielectric slab 
shown in Fig. 7b. In that case, E has only one component: E v = E(x), 
H z = (l/iuno)dE/dx, and H x — —(kJu>n )E. Taking into account 
Sf = E y and H} = H z (see Ref. 1), we obtain 

Ak z = l(EdE p /dx) - (E p dE/dx)l/(2k z f +t °E 2 dx\ , (55) 

where we have assumed that E, differs significantly from E only near 
the perturbing slab. The unperturbed field is, for — D m < x < 0, 

E = exp (— kx), (56) 

and the perturbed field is that given in (45) 

E p = exp ( — kx) + r exp {kx), (57) 

where 

k - {k\ - k 2 ,)K 

The amplitudes in (56) and (57) are so chosen that E p ttE for large 
negative x, e.g., x = —D m . 
We first evaluate 

(EdE p /dx) - (E p dE/dx) = 2r/c, (58) 

where we have used (56) and (57). Note that the result (58) is inde- 
pendent of x (for —D m < x < 0). Substituting (58) in (55), the 
imaginary part of k, is found 

k ai = r iK Rexp (-2 K D m ), (59) 

where r, denotes the imaginary part of r, evaluated in Appendix B. 
We have introduced in (59) the field strength parameter 

R = (k z j +X EHxX\ (60) 

TRANSVERSE COUPLING IN FIBER OPTICS IV 1449 



In the above definition of R, the field is assumed to be unity at the 
guide-cladding boundary. For a dielectric slab, the value of R is given 
in (12). Equation (59) is the result used in the main text. 
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